Charmed-Meson Decay Constants in Three-Flavor Lattice QCD 
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We present the first lattice QCD calculation with realistic sea quark content of the Z) + -meson 
decay constant f D + . We use the MILC Collaboration's publicly available ensembles of lattice gauge 
fields, which have a quark sea with two flavors (up and down) much lighter than a third (strange). 
We obtain f D + = 201±3±17 MeV, where the errors are statistical and a combination of systematic 
errors. We also obtain fn B = 249 ± 3 ± 16 MeV for the D s meson. 

PACS numbers: 13.20.Fc,12.38.Gc 



Flavor physics currently plays a central role in elemen- 
tary particle physics Q. To aid the experimental search 
for physics beyond the standard model, several hadronic 
matrix elements must be calculated nonperturbatively 
from quantum chromodynamics (QCD). One of the most 
important of these is the decay constant of the B meson 
/b Hi- Any framework for calculating fs should, there- 
fore, be subjected to stringent tests, and such a test is a 
key aim of this Letter. 

The most promising method for these nonperturbative 
calculations is numerical lattice QCD. For many years 
the results suffered from an unrealistic treatment of the 
effects of sea quarks. In the last few years, however, this 
obstacle seems to have been removed: with three flavors 
of sea quarks lattice QCD now agrees with experiment for 
a wide variety of hadronic quantities [j| . This validation 
of lattice QCD has been realized, so far, only for so-called 
"gold-plated" quantities: masses and matrix elements of 
the simplest hadronic states. Note, however, that many 
of the hadronic matrix elements relevant to flavor physics 
are in this class, including fg. 

The challenges in computing fg are essentially the 
same for the D + -meson decay constant Jd+- Experi- 
ments have observed the leptonic decay D + — ► l + vi, but 
not B + — ► l + vi. One can, thus, determine |V c d|/£i+, 
where V c d is an element of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix. Taking | V^d | from elsewhere, 
one gets fo+- In 2004 the CLEO-c Collaboration mea- 



sured /£>+ with a 20% error jj] , and a more precise mea- 
surement is expected soon. 

This Letter reports the first lattice-QCD calculation of 
fu+ with three flavors of sea quarks We find 



Id* 



201 ± 3 ± 6 ± 9 ± 13 MeV, 



(1) 



where the uncertainties are statistical, and a sequence of 
systematic effects, discussed below. We also obtain the 
decay constant of the D s meson, 



f Ds = 249 ± 3 ± 7 ± 11 ± 10 MeV. 



(2) 



The second result is more precise than a recent lattice- 
QCD calculation with the same sea quark content but 
non-relativistic heavy quarks, which found fu s = 290 ± 
20 ± 41 MeV 0. These results are more reliable than 
older calculations Q because we now incorporate (three) 
sea quarks and, for fu+, also because the light valence 
quark masses are smaller than before. 

These results test the methods of Ref. Q because they 
are predictions. The input parameters have been fixed 
previously 0, El E EI , an d, once comparably precise 
experimental measurements become available, one can 
see how Eqs. Q and @ fare. Indeed, this work is part of 
a program to calculate matrix elements for leptonic and 
semileptonic decays hfl H2I Il3j | , neutral-meson mixing, 
and quarkonium [111 So far, these lattice QCD cal- 
culations agree with experiment for the normalization of 
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TABLE I: Notation for quark masses used in this Letter. 



m 


Description 




Remark 


m c 


Charmed quark 


From 


m Ds [10, nj 


in s 


Physical strange quark 


From 


m\ [8] 


m u 


Physical up quark 


m u = 


m s /45.5 [9] 


m d 


Physical down quark 


m d = 


m s /19.6 [9J 


rrih 


Simulation's heavier sea quark 


m h » 


l.lm s 


mi 


Simulation's lighter sea quark 


0.1m s 


< m; < 0.8m s 


m q 


Simulation's light valence quark 


0.1m s < rriq < m s 



Z?-meson semileptonic form factors 0,^3,^3- They also 
have predicted correctly the form- factor shape |12t 1 1 7| , 
as well as the mass of the B c meson 0, 0] . 

In this set of calculations we use ensembles of un- 
quenched latti ce g auge fields generated by the MILC Col- 
laboration |9|,ll9|) with lattice spacing a — 0.175, 0.121, 
and 0.086 fm. The key feature of these ensembles is that 
they incorporate three flavors of sea quarks, one whose 
mass is close to that of the strange quark, and two with 
a common mass taken as light as possible. 

For the sea quark and light valence quark we use the 
"Asqtad" staggered- fcrmion action |20j. Several differ- 
ent quark masses appear in this calculation; for conve- 
nience, they are defined in TablcQ] At a — 0.175, 0.121, 
and 0.086 fm there are, respectively, 4, 5, and 2 ensem- 
bles with various sea quark masses (m^m.^) 0,0- The 
larger simulation mass, rrih is close to the physical strange 
quark mass m s . The light pair's mass m; is not as small 
as those of the up and down quark in Nature, but the 
range 0.1m s < mi < 0.8m s suffices to control the ex- 
trapolation in quark mass with chiral perturbation the- 
ory (%PT). For carrying out the chiral extrapolation, it 
is useful to allow the valence mass m q to vary separately 
from the sea mass [2lj]. At a = 0.175, 0.121, and 0.086 fm 
we have, respectively, 6, 12, and 8 or 5 values of the va- 
lence mass, in the range 0.1m s < m q < m s . 

A drawback of staggered fermions is that they come 
in four species, called tastes. The steps taken to elim- 
inate three extra tastes per flavor are not (yet) proven, 
although there are several signs that they are valid. Cal- 
culations of f D + and fo s are sensitive to these steps: 
if Eqs. and © agree with precise measurements, it 
should be more plausible that the techniques used to re- 
duce four tastes to one are correct. 

For the charmed quark we use the Fermilab action for 
heavy quarks [j^]. Discretization effects are entangled 
with the heavy-quark expansion, so we use heavy-quark 
effective theory (HQET) as a theory of cutoff effects [23j . 
This provides good control, as discussed in Ref. 0, and 
the framework has been tested with the (successful) pre- 
diction of the B c meson mass 0]. Nevertheless, heavy- 
quark discretization effects are the largest source of sys- 
tematic error in /jj a , and the second-largest in fjj+. 

The decay constant /n , for a D q meson with light 



valence quark q and momentum p^ , is defined by (2f 



(3) 



where = 97^ 75 c is an electroweak axial vector cur- 

t> q = fD q y/mr^ emerges di- 



rent. The combination 

rectly from the lattice Monte Carlo calculations. As 
usual in lattice gauge theory, we compute two-point cor- 
relation functions C 2 (t) = (O]j g (t)On t (0)), C A {t) = 
(Ai{t)0 D q {fy) 1 where Oo q is an operator with the quan- 
tum numbers of the charmed pseudoscalar meson, and 
Ai is the (lattice) axial vector current. The operators are 
built from the heavy-quark and staggered-quark fields as 
in Ref. (2(|. We extract the D q mass and the ampli- 
tudes (D\O Dq \0) and (0\A 4 \D) from fits to the known t 
dependence. Statistical errors are determined with the 
bootstrap method, which allows us to keep track of cor- 
relations. 

The lattice axial vector current must be multiplied by 
a renormalization factor Zg* . We write |27| Z A ^i = 

p A cq {Zy^Z v <ii ) 1 ^ 2 , because the flavor-conserving renor- 
malization factors Zycc and Zym are easy to compute 
nonperturbatively. The remaining factor p A ^i should be 
close to unity because the radiative corrections mostly 
cancel (2^. A one- loop calculation gives [2^| p A <=i = 
1.052, 1.044, and 1.032 at a = 0.175, 0.121, and 0.086 fm. 
We estimate the uncertainty of higher-order corrections 
to be 2a s {pA a 4 q — 1) ~ 1.3%; a s is the strong coupling. 

The heart of our analysis is the chiral extrapolation, 
from the simulated to the physical quark masses. It is 
necessary, and non-trivial, because the cloud of "pions" 
surrounding the simulated D q mesons is not the same 
as for real pions. With staggered quarks the (squared) 
pseudoscalar meson masses are 



{m a + m b )n ■ 



a 2 A 



(4) 



where m a and nib are quark masses, p is a parameter 
of xPT, and the representation of the meson under the 
taste symmetry group is labeled by £ = P, A, T, V, I [30| . 
A symmetry as m Q ,m(, — > ensures that Ap = 0. The 
"pion" cloud in the simulation includes all these pseu- 
doscalars. 

According to next-to-leading order %PT the decay con- 
stant takes the form 

4> q = $ [1 + Af q (m qi mi,m h ) + p q (m g ,mi,m h )] , (5) 

where $ is a quark-mass-independent parameter. Af q 
arises from loop processes involving light pseudoscalar 
mesons, and p q is an analytic function. To obtain them 
one must take into account the flavor-taste symmetry of 
the simulation 30] and the inequality (in general) of the 
valence and sea quark masses [Hi. One finds [HI 



Af q 



l + 3g 2 



[h q + h q + a 2 (S' A hf + S' v h^)], (6) 
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where f n f» 131 MeV is the pion decay constant, g is the 
D-D* -it coupling [32, and S' A , S' v parametrize effects that 
arise only at non-zero lattice spacing j^B]. The terms h q , 
h q , hg, and h^ are functions of the pseudoscalar meson 
masses. The last two, h q and hX , are too cumbersome 
to write out here. It is instructive to show the other two, 
h q and h ! q , when m q = mi or m^: 

K = ^E«»« [^(M^)+I(M^)] , (7) 
h\ = -i/(M2 ;J ) + i7(M n 2 >7 ), (8) 
^ = —I(M% hI ) + |7(M^j), (9) 

where 7(Af 2 ) = M 2 In Af 2 /A 2 (with A x the chiral scale), 
and M 2 j = (Mfij + 2M% hjI )/3. The term h\ receives 
contributions only from taste-singlet mesons (represen- 
tation /). The term h q receives contributions from all 
representations, with multiplicity = 1,4,6,4,1 for 
£ = P, A, T, V, I, respectively. The analytic function is 



p q = (2rrn + m h )h{A x ) + m q fc{A x ) + 0(a 2 ), (10) 

where f\ and f% are quark-mass-independent parameters. 
They are essentially couplings of the chiral Lagrangian, 
and their A x dependence must cancel that of Af q . This 
specifies 0(a 2 ) terms proportional to /i and f%, which 
can be removed after our fit. We estimate the remaining 
0(a 2 ) effects of light quarks to be small: around 4% at 
a = 0.121 fm and 1.4% at a = 0.086 fm. 

The salient feature [3^ of the chiral extrapolation of (fi q 
is that Afq contains a "chiral log" I(2m q fi) ~ m q \nm q , 
which has a characteristic curvature as m q — > 0. Equa- 
tions I©-© show that the chiral log is diluted by dis- 
cretization effects, because a 2 A^ ^ for £ ^ P. 

We can now discuss how we carry out the chiral ex- 
trapolation. Recall that we compute <j) q for many com- 
binations of the valence and light sea quark masses. At 
each lattice spacing, we fit all results for (fi q to the mass 
dependence prescribed by Eqs. ijlJl- ljTUI) . Of the twelve 
parameters, eight — ji, the four non-zero A^, f„, S' A , and 
S' v — appear in the x?T for light pseudoscalar mesons. 
We constrain them with prior distributions whose central 
value and width are taken from the x?T analysis of pseu- 
doscalar meson masses and decay constants on the same 
ensembles of lattice gauge fields Q. The rest — $, g 2 , /i, 
and fi — appear only for charmed mesons. We constrain 
g 2 to its experimentally measured value, within its mea- 
sured uncertainty jl^l- Thus, only three parameters — 
fx, and /2 — are determined solely by the <fi q fit. To obtain 
physical results we reconstitute the fit setting the light 
sea quark mass mi — + (m u + md)/2, and = S' A v — 0. 
For <fid {<fis) we se t the light valence mass m q md (m s ). 

To isolate the uncertainties of the chiral extrapolation 
from other sources of uncertainty, we consider the ratio 
Rq/s — (fiq/^s- Figuren snows Rq/s at o = 0.121 fm as a 
function of m q /m s , projected onto m q = mi. The gray 



0.90 r 




FIG. 1: Chiral extrapolation of R q / S at a — 0.121 fm. Data 
points show only statistical errors, but the systematic error 
of fitting is shown at left. 




FIG. 2: Dependence of <fis on a 2 . Circles result from removing 
the 0(a 2 ) pieces in Eq. 1101 : squares omit this step. 



(red) curve is the result of the full fit of <fi q to the separate 
sea- and valence-mass dependence. The black curve, and 
the extrapolated value at m q /m s = 0.05, results from 
setting A ^ = S' A v — when reconstituting the fit. At 
the other lattice spacings we obtain similar results. 

The precision after the chiral extrapolation is, however, 
a bit illusory. We tried several variations in the fit proce- 
dure: fitting the ratio directly; adding terms quadratic in 
the quark masses to Eq. ifTfljl : variations in the widths of 
the prior constraints of the parameters. When these pos- 
sibilities are taken into account, the extrapolated value 
of Rd/s varies by 5%, which we take as a systematic un- 
certainty. This variation could be reduced with higher 
statistics at the lightest sea quark masses. 

The lattice spacing dependence of <fi s — /d sA /tod7 is 
shown in Fig- El The (blue) circles are the main results. 
In a preliminary report of this work &J the 0(a 2 ) terms 
in <fi s were not removed. The (red) squares illustrate the 
effect of omitting this step. As one can see, the effect is 
small at a = 0.086 fm, but it is the main reason whythe 
results in Eqs. and (J2J are smaller than in Ref. |5j. 

The %PT expressions for (fi q assume that the D q me- 
son is static. Since its mass is around 1900 MeV and 
the pseudoscalars are a few hundred MeV, this is a good 
starting point. Some corrections to this approximation 
can be absorbed into the fit parameters, with no real 
change in the analysis. A more interesting change arises 
in the one-loop self-energy diagrams, for which the func- 
tion J(M 2 ) is modified, and depends on mu* — mu as well 
as M. By replacing our standard extrapolation by one 
using the modified function, we estimate the associated 
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TABLE II: Error budget (in per cent) for R d / 


j, <f> s , <; 




source 


Rd/s 


4>s 




statistics 


0.5 


1.4 


1.5 


input parameters a and m c 


0.6 


2.8 


2.9 


higher-order p A ^g 





1.3 


1.3 


heavy-quark discretization 


0.5 


4.2 


4.2 


light-quark discretization and xPT fits 


5.0 


3.9 


6.3 


static xPT 


1.4 


0.5 


1.5 


finite volume 


1.4 


0.5 


1.5 


total systematic 


5.4 


6.5 


8.5 



error to be 1.5% or less. Finite-volume effects also mod- 
ify I(M 2 ): based on our experience with / w and fx 
and on continuum xPT pjjj ]. we estimate a further error 
of 1.5% or less. 

Although xPT is able to remove (most of) the light- 
quark discretization errors, heavy-quark discretization ef- 
fects remain. We estimate this uncertainty using HQET 
as a theory of cutoff effects j^, 01 • P° arrive at a nu- 
merical estimate, one must choose a typical scale A for 
the soft interactions; we choose A w 500-700 MeV. We 
then estimate a discretization uncertainty of 2.7-4.2% at 
a = 0.086 fm. Similarly, the results at a = 0.121 fm are 
expected to lie within 1-2% of those at a — 0.086 fm. 

Because we cannot disentangle heavy- and light-quark 
discretization effects, to quote final results we average the 
results at a — 0.086 and 0.121 fm. We then find 

R d/S = 0.786(04) (05) (04) (42) (11) 
(j> s = 0.349(05) (10) (15) (14) GeV 3/2 , (12) 

which are the principal results of this work. The uncer- 
tainties (in parentheses) are, respectively, from statistics, 
input parameters a and m c , heavy-quark discretization 
effects, and chiral extrapolation. A full error budget is 
in Table ILTl all uncertainties are reducible in future work. 
The results for f D + and /jj s in Eqs. (JTJ and J5J are ob- 
tained via f Ds = (j> s /^/m Ds , f D + = Rd/s4's/y / m D +, by 
inserting the physical meson masses. 

Present experimental measurements, fu+ = 202±41± 
17 MeV P|, f Ds = 267 ±33 MeV jH, are not yet precise 
enough to put our results in Eqs. JU and J5J to a strin- 
gent test. The anticipated measurements of fo+ and, 
later, fu t from CLEO-c are therefore of great interest. 
If validated, our calculation of fjj+ has important impli- 
cations for flavor physics. For B physics it is crucial to 
compute the decay constant fg. To do so, we must sim- 
ply change the heavy quark mass. In fact, heavy-quark 
discretization effects, with the Fermilab method, are ex- 
pected to be smaller, about half as big. 

We thank the U.S. National Science Foundation, the 
Office of Science of the U.S. Department of Energy, Fer- 
milab, and Indiana University for support, particularly 



for the computing needed for the project. Fermilab is 
operated by Universities Research Association Inc., un- 
der contract with the U.S. Department of Energy. 

Note added: After this Letter was submitted, the 
CLEO-c Collaboration announced a new measurement, 
f D + = 223 ± 16± 7 g MeV 
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